UHECR Production by a Compact Black-Hole Dynamo: Application to Sgr A* 
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The possibility that the excess cosmic ray (CR) flux near 10 18 eV, reported recently by the 
AGASA group, is due to a compact black hole dynamo associated with the Sgr A* source is consid- 
ered. Under the assumption that the Galactic center black hole rotates with nearly maximal spin, 
and that the magnetic field threading the horizon is in rough equipartition with matter accreted by 
the hole, the spectra and total fluxes of accelerated CRs and their associated curvature emission 
depend only on accretion rate. For the accretion rate estimated on the basis of observations of stellar 
winds near Sgr A* (~ 1CT 3 Eddington), the maximum proton energy achievable by this mechanism 
is of the order of 10 18 eV, and the corresponding CR power is ~ 10 41 qca erg s _1 , where acR « 1 
is the CR production efficiency. The corresponding spectrum of curvature photons peaks at around 
a hundred GeV, with a total luminosity comparable to the CR power released. For much lower 
accretion rates both the CR and gamma-ray fluxes predicted are completely insignificant. Upcom- 
ing gamma-ray experiments, such as CANGAROO, HESS and GLAST, can be used to probe the 
parameters of this system. 
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I. INTRODUCTION 



Recent analysis of data from AGASA fl|] and SUGAR |2| have indicated an excess of cosmic ray intensity over a 
narrow energy range around 10 18 eV, that has been interpreted as due to a strong CR source in the direction of the 
Galactic center. The analysis by the AGASA group suggested a rather extended source in the general vicinity of the 
Galactic center, with the most significant excess obtained for a signal beam size of about 20°. Close examination of 
the SUGAR data that was motivated by the results of the AGASA group, confirmed the existence of an excess CR 
flux near the direction of the Galactic center, but indicated a signal that appears to be consistent with a point source 
offset by about 7.5° from the true Galactic center. Clay ||| suggested that the pointlikc excess is caused by neutrons 
produced as a result of conversions of particles in a target located outside their acceleration site; in order to reach 
Earth before decaying, the neutrons should have energies in excess of ~ 10 18 eV. This naturally explains the absence 
of an excess at energies below 10 18 eV. The absence of an excess at energies well above that is attributed to an upper 
limit on the acceleration energy in the source. Clay (see also ref. also proposed that the propagation of the charged 
(non-converted) particles through the Galactic magnetic field might give rise to the more extended source that appears 
in the AGASA data. Using detailed numerical simulations and reasonable models for the Galactic magnetic field he 
studied the propagation of protons from the vicinity of the Galactic center to Earth, and concluded that a source 
of protons near the Galactic center would be unobservable at energies below 10 18 eV, but should produce a halo of 
energy dependent size and location above this energy. The details depend on the magnetic field model adopted, but 
in general there should be a systematic shift of the source image to the north, owing to the regular field component. 

The above considerations motivated us to explore the application of a model, developed recently to explain the origin 
of the highest energy cosmic rays detected ||, to the Galactic center black hole, Sgr A*. In this model the ultra-high 
energy cosmic rays (UHECRs) are accelerated by the electric potential difference generated by spinning supermassive 
black holes associated with dormant AGNs. The model assumes that those black holes were either formed with critical 
angular momentum or spun up during an earlier phase when the AGNs were active, and are presently rotating with 
nearly maximal spins. The basic picture envisaged is that a small number of particles injected into the gap are being 
accelerated to ultra-high energies during episodes when the emf induced by the hole dynamo is not shorted out; in 
this picture the rotational energy of the hole is liberated essentially in the form of UHE particles and their associated 
curvature emission ||, rather than powerful jets, as seen in blazars. This is consistent with the expectation that 
spontaneous breakdown of the vacuum should not occur in those systems Q . For the range of parameters relevant to 
UHECR production, the curvature losses suffered by the accelerating particles are severe, suppressing the maximum 
energy they can achieve substantially [^[0]. The curvature photons are emitted predominantly in the TeV band, with 
an average flux per UHECR source that exceeds the detection limit of current TeV experiments || , and is conceivably 
highly variable. The total power needed to be extracted from an average UHECR source in order to account for the 
measured CR flux above the GZK cutoff is rather small, less than about 0.1% of the maximized Blandford-Znajek 
(BZ) power. Why this process is so inefficient as compared with the process of jet formation in blazars is yet an open 
issue. As argued by Punsley & Coroniti Q the BZ power extracted from a BH should be governed by the process of 
plasma injection in the ergosphere. So perhaps in the UHECR sources, for which spontaneous vacuum breakdown as 
in blazars is not expected, plasma injection is considerably suppressed, resulting in a much lower efficiency. 

Recent numerical simulations Q suggest that the accretion process and magnetic field structure in the vicinity 
of the horizon should be non stationary, owing to rapid magnetic field reconnection. This would probably lead to 
appreciable complications of the (over) simplified model adopted in ref [j5| . For instance the location of the gap and, 
perhaps, the voltage drop across it might change with time, as well as the injection of seed particles into the gap. 
How should this affect the picture described above is unclear at present. Despite of all the complications anticipated, 
we shall proceed by using the simple model j^] outlined above. 

II. APPLICATION TO SGR A* 

There are essentially three important parameters that determine the CR and gamma-ray fluxes produced by the 
black hole dynamo: the mass of the hole, its angular momentum, and the strength of the magnetic field threading the 
horizon. In the case of Sgr A* the mass is known to a reasonably good accuracy from dynamical measurements. The 
other two parameters, however, arc highly uncertain. Upper limits on the power and spectral peak energy of emitted 
cosmic rays can be obtained if the hole angular momentum is taken to be near its critical value, and the magnetic field 
is assumed to be in equipartition with the matter accreted into the black hole. Under these circumstances the results 
depend only on one parameter, namely the accretion rate m, henceforth measured in Eddington units (LEdd/c 2 )- 
Below we adopt the estimate of Eckart & Genzel |1Q] for the black hole mass: M — 2.5 x 10 6 M©. 
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The rate at which the Galactic black hole is accreting is unfortunately unknown. The total luminosity associated 
with Sgr A* is of order 10 37 erg s _1 (see ref. |ll]]). If we naively relate this luminosity to accretion from a standard 
thin disk with say 10 % efficiency, we obtain to ~ 10 -7 . This value is well below those estimated by Melia JO] 
(to ~ 1Q~ 15 ) and Genzel et al. Jl3| (m ~ 10~ 3 ), based on observations of stellar winds in the vicinity of Sgr A* and 
assuming Bondi accretion (that might be irrelevant if the infalling matter possesses angular momentum). Narayan 
et al. constructed an ADAF model for Sgr A* in an attempt to reconcile the low luminosity of Sgr A* with the 
relatively high accretion rates estimated by Melia and Genzel et al. A viable fit of the spectrum of Sgr A* by the 
ADAF model (but cf. ref. jl4| for a different interpretation) appears to be consistent with the measured mass of 
Sgr A*, with an accretion rate of to ~ 10~ 3 , and with an equipartition magnetic field, provided the parameter S is 
sufficiently small (5 ~ 10 -3 ), where 5 represents the fraction of turbulent energy that is tapped for electron heating. 
However, it has been subsequently shown Jl5j that models with large S and considerable mass loss rate via winds [ |l6| , 
for which only a small fraction of the inflowing matter ultimately reaches the black hole, are also in agreement with 
observations. Thus, it seems that there is, at present, a large uncertainty in the value of to, but that it is likely to lie 
in the range 10~ 7 < to < 10 -3 . 

The strength of the equipartition magnetic field depends, quite generally on the details of the inflow pattern in 
the vicinity of the horizon. A rough estimate gives B4 ~ Ql(rn/M§) 1 l' 2 , where Me = A//1O 6 M ||. This is in good 
agreement with the mean equipartition field calculated using the self-similar ADAF model by Narayan and Yi [fj"7| . 
For Sgr A* (M 6 = 2.5) wc obtain 

B 4 ~ 38to 1/2 . (1) 

The emf generated by a rotating black hole is given approximately by AV ~ 4.5 x 10 17 (a/M)(h/ ' R g ) 2 B 4 MQ V. 
Taking a — M, h — R g , M 6 — 2.5, and using eq. (Q) yields 

AV ~5x 10 19 to 1/2 V. (2) 

For the range of to considered here (10 -7 to 10~ 3 ) the emf given by the last equation lies in the range between about 
10 16 and 1.5 x 10 18 Volt. 

The energy change per unit length of an accelerating particle of charge Z and energy e = jrriiC 2 can be expressed 
as §, 

de _ eZAV 2 e 2 Z 2 7 4 
ds h 3 p 2 

where p is the mean curvature radius of magnetic field lines in the gap. The first term on the r.h.s of the last equation 
accounts for the energy gain owing to the acceleration by the gap electric field, and the second term accounts for the 
losses due to curvature radiation. Under the assumption that p is independent of the particle's energy, eq. (|J) can be 
solved analytically. The maximum energy achievable, e max , can then be obtained by integrating eq. (^) over the gap 
(i.e., from s = to s = h), and is given implicitly by, 

lniif + 2 tan -1 (a) = 477 -1 . (4) 
1 — x 

Here 77 = 3p,(a/M)- 3 / 4 M~ 1/2 Z- 5 /' i B~[ 3/4 (p/R g ) 1 / 2 (h/R g )~ 7 /' 1 , with p, being the mass of the ion in units of the 
proton mass, is the suppression factor defined in eq. (5) of ref. Q for a — M, and x = e max /e max , where e m ax, 
given explicitly in eq. (4) of ref. ||, is the maximum acceleration energy in the limit of large suppression (i.e., strong 
curvature losses), as can be readily seen from eq. (Q) wherefore x ~ 1 when 77 << 1. In the opposite limit, ?y >> 1, 
the solution to eq. (||) is given to a good approximation by e max ~ eZAV, as expected. For the parameters adopted 
above we find 

n ^ o.i2pZ- 5/ Sh~ 3/8 , (5) 

where it has been assumed that p = h = R g . 

The maximum rotational power that can be extracted from a Kerr black hole is given roughly by Lbh ~ 
10 40 (a/M) 2 M|_B| erg s -1 . We suppose that a fraction atcR « 1 of this power is liberated in the form of CRs. 
Taking again a — M, and using eq. (|l|), we find that the CR power released by the Galactic BH can be written as 

L C r ~ 10 44 a C i?m ergs -1 . (6) 
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The relative production efficiency of curvature photons, defined as the ratio of the total curvature loss per nucleus 
and the maximum energy gain, can now be expressed as 

k = * 0.5Z 5 5| M - 4 (e max /e^AF) 3 ~ 8 x 10 4 Z 5 ^m^ie^/eZAV) 3 , (7) 

^max 

where the dimensionless parameter, e max /eZAV, can be computed by employing eqs. ^ and (||) once m is specified. 
The limit of large suppression corresponds to k >> 1, whereas k « 1 corresponds to small curvature losses. To a 
good approximation then the luminosity associated with curvature emission is given by L 7 ~ kLcr- The curvature 
spectrum will peak at an energy || 

e 7 - 60{Z/nfBl = 3 x 10 6 (Z//^) 3 m 3/2 GeV. (8) 

We now use the results derived above to estimate the luminosity and spectral peak energy of emitted CRs and 
curvature photons, for different accretion rates corresponding to different models of Sgr A*. Consider first the 
possibility of standard accretion with a very low rate. Adopting rh = 1CP 7 (see above), eq. (pi) yields rj ~ 51 for 
protons and r\ ~ 47 for iron. Substituting the latter result into eq. (Q), one finds e max — 1.5 x 10 eV for protons and 
e max ~4x 10 17 eV for iron. The spectrum of curvature photons should peak at around 100 keV if the accelerated 
particles are predominantly protons and 10 KeV if iron. From eq. (0) we obtain a CR power of Lcr — 10 37 ac_R erg 
s _1 , and from eq. (0) we obtain a radiative efficiency of k ~ 10 -5 ' 5 for both protons and iron and a corresponding 
gamma-ray luminosity of L 7 ~ 3 x 1Q 31 acR erg s _1 . We therefore conclude that for such a low accretion rate the 
fluxes of CRs and curvature photons produced by the dynamo mechanism are completely negligible. 

We consider next the situation whereby the radio - through X-ray spectrum of Sgr A* is produced by an ADAF. 
Adopting the accretion rate obtained from the standard ADAF model |0|, viz., m ~ 10~ 3 , and repeating the above 
calculations we find: e max ~ 10 18 eV for protons and e max — 2.5 x 10 19 eV for iron, a CR power of Lcr — 10 41 acfl 
erg s _1 , a radiative efficiency of k ~ 1 for both protons and iron with a corresponding gamma-ray luminosity 
Lj ~ Lcr, and peak energy for the curvature spectrum of approximately 100 GeV or 10 GeV, depending on whether 
the accelerated particles are, respectively, protons or iron. We stress that the estimate of the peak energy is highly 
uncertain because of the strong dependence of e 7 on the strength of the magnetic field (see eq. Q). Taking for 
illustration acR = 10~ 3 , which is roughly the value inferred for the dormant AGNs (|||, yields Lcr ~ L 1 ~ 10 38 erg 
s _1 . The corresponding gamma-ray flux at Earth (S 7 ~ 10~ 8 erg cm -2 s _1 ) is several orders of magnitude above the 
threshold sensitivity anticipated for the HESS pl| and CANGAROO [jl9| imaging atmospheric Cerenkov telescopes 
in the southern hemisphere. The photon detections to be obtained with GLAST during a normal one-year scan mode 
are expected to be on the order of a thousand counts or more in this case (Thompson 2000, personal communication). 

In addition to the energy losses corresponding to curvature radiation we must also address the possible impact of 
inelastic pair and pion producing collisions with the environmental electromagnetic radiation near the hole (see ref. [Q). 
We estimate that, during the acceleration of a proton to 10 18 eV by the putative dynamo associated with Sgr A*, there 
should be negligible drag arising from inelastic pair producing collisions with ambient photons. In particular, since the 
Lorentz gamma factor is never larger than 10 9 during the acceleration process, only submillimeter (wavelength < 1.2 
mm) ambient photons are involved in electron-positron pair production. From the estimated SED of Sgr A* [pTpil we 
conclude that the corresponding number density of suitable submillimeter photons implies a radiation length about 
two orders of magnitude larger than the intrinsic size of the radio source [[22| . Finally, it is important to note that 
collisions involving pion generation are of substantially higher inelasticity than those involving pair production |2^j 
and thereby potentially much more catastrophic for the putative dynamo action. However, for inelastic collisions with 
photo-pion production, only those ambient photons of wavelength shorter than 8 microns are involved; the associated 
collision mean-free-path in this instance is then estimated to be even larger than the radiation length characterizing 
pair production. 



III. DISCUSSION 



In this paper we have considered the application of a compact black hole dynamo mechanism to CR production by 
Sgr A*. We have argued that if the rate at which matter is accreted by the black hole is on the order of that estimated 
on the basis of observations of stellar winds near Sgr A*, then protons of maximum energy ~ 10 18 eV and total power 
~ 10 38 erg s^ 1 (assuming CR production efficiency of 10~ 3 ) can be produced by the dynamo mechanism, provided 
the black hole angular momentum is near its critical value, and assuming equipartition between the magnetic field 
threading the horizon and inflowing matter. An associated gamma-ray emission arising from curvature acceleration, 
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with a luminosity comparable to the total CR power released, and with a spectrum that peaks sharply at around 
100 GeV is predicted under these conditions, and motivates observations with future missions. Positive detection of 
curvature photons will have important implications for the process of accretion into the black hole. 

Whether the above scenario can explain the excess CR flux reported by the AGASA group recently, even in the case 
of large accretion rate, remains to be checked; at these energies propagation effects appear to be quite sensitive to the 
Galactic magnetic field model adopted || and should be accounted for properly. The energy flux of the apparent point 
source offset from the galactic center observed with SUGAR is 3 x 10~ n ergs s _1 cm~ 2 B. Assuming its distance 
from Earth to be comparable to that of Sgr A*, yields a luminosity of ~ 2 x 10 35 ergs/(s en?). If the SUGAR source 
is associated with energetic neutrons produced in a target there by the conversion of UHECR protons (or Fe nuclei) 
originating from Sgr A*, then given the error circle of SUGAR, it implies that the target size should lie in the range 
between 30 to 300 pc. 

Finally, it is interesting to note that the galactic center region is in fact a known gamma-ray source |24|, observed 
over the band 30 MeV to 10 GeV with the EGRET instrument on the Compton GRO to exhibit a hard power law 
spectrum that steepens above 2 GeV and where the inferred energy flux within the total measured band is 3 x 10~ 9 
ergs s _1 cm~ 2 ; however, the emission volume is likely extended, and the portion arising from a possibly major point 
source is still uncertain (an issue that should be much better addressed with the improved angular resolution to be 
provided by GLAST). If we naively associate the break with the peak energy of curvature photons, we obtain from eq. 
(||) rh ~ 10~ 45 . From eqs. (||) and (0) we obtain a luminosity of L 7 ~ 10 3S acR erg s _1 . Assuming that about half 
the flux is emitted by the point source 0], we obtain a CR production efficiency of acR < 0.1, which seems rather 
high compared with the efficiency inferred for the dormant AGNs. However, since e 7 oc 5 3 , as seen from eq. (||), 
even relatively small deviations from equipartition would lead to considerably lower peak energy and, consequently, 
higher values of rh and lower values of acR- Alternatively, the potential curvature radiation from Sgr A* might 
well be a separate gamma ray component from that detectable with EGRET (e.g., at a much higher energy, such as 
characteristic of the case illustrated here) . 

We thank Michael Loewenstein, Gustavo Medina- Tanko and David Thompson for enlightening discussions. 
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